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Abstract: We show that quantum particles constrained to move along curves undergoing cyclic 
deformations acquire, in general, geometric phases. We treat explicitly an example, involving particular 
s21 ' deformations of a circle, and ponder on potential applications. 
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1 Introduction 

You leave a wire on a table in a locked room, and rumours soon have it that somebody entered the room, 
played with the wire, and left it at the exact same initial position — is it possible to detect the intrusion? 

Geometrical phases may appear in quantum systems the hamiltonian H of which depends on a set 
of external parameters ^'^, H = H{^). After a cyclic change of the parameters, along a loop C in 
parameter space, slow enough to guarantee that the system remains in an instantaneous eigenstate of H, 
the wavefunction acquires a phase that, apart from the standard dynamical part of the "energy times 
time" type, contains a contribution that only depends on C — hence the term "geometric" 1 . Observable 
effects are obtained if the system starts out in a suitable superposition of energy eigenstates, each of which 
acquires a different geometric phase — numerous experiments have confirmed these theoretical predictions 
(see, e.g., [HI El [HI [HI EH [l6| ) . 

On a different vein, the dynamics of quantum particles constrained to move along a curve in three 
dimensional euclidean space, by a steep confining potential in the plane normal to the curve, is governed 
by an effective hamiltonian that depends on the curvature k(s) and the torsion t{s) of the curve, where 
s denotes arclength [71 [H [T71 [TU] (see also [HI [HI [3 [HI E] ) • Thus, the Fourier coefficients, for example, of 
K and r, can be considered as external parameters of the effective hamiltonian. Slow, cyclic changes in 
the shape of the curve may then lead to the appearance of geometric phases, and it is the confirmation 
of this possibility that we report on here. Thus, the opening question is answered to the affirmative, at 
least for a playful enough intruder (see below). Apart from limiting what intruders can get away with, 
the effect pointed out here opens up a wide arena for experimentation, both gedanken and real — we 
ponder on some possibilities in the concluding section. 

2 Geometric phases 

Consider a hamiltonian H^, as above, where the ^'s are varying with time, tracing out a loop C in ^-space. 
Suppose the physical system governed by H^ starts out, at time i = 0, in the nondegenerate eigenstate 
\n, ^o) 7 which satisfies H^^ \n, ^o) = E^" {n, ^o) • Then, in the adiabatic approximation, in which the change 
in H is slow (in the time scale set by the energy difference of neighboring eigenstates), the system's state 
at time t is 

|n,t)-e— "(*)+^^"W|n,6), (1) 

where \n,£,t) is an instantaneous eigenket of the hamiltonian, H^^\n,S,t) — E^\n,^t)i the phase a„(i) = 
Jp dr E^ is the expected dynamical one, and 

is the geometrical phase — the latter form shows that it is time reparametrization invariant (V^ denotes 
the gradient in ^-space). For a loop C in ^-space, Stokes' theorem may be used to cast ^ in the form 

7n(C)= I d^d^K'^AL (3) 

Js 

where 

K^^l^7{dA{n,mdB\n,£.)) (4) 

is the Berry curvature, and S is any two-dimensional patch with C as its boundary [I] 1141 [5] (see 
also 0[H1)- 

When the initial state \n,t = 0) is d-fold degenerate, the geometric phase factor e'^"*^'^ generalizes to 
a unitary matrix [20] (the Wilczek-Zee effect) 

Unit) = Pe^pliJ' d^Ai\ , (5) 
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where P exp is a path-ordered exponential, and A^^ is a hermitean d by d matrix with entries 

{Ai)ab^i{n,b;^\V^\n,a;0 (6) 

(a, b = 1, . . . ,d, range over the degenerate subspace). 

3 Particles Constrained on Curves 

Consider a quantum particle constrained to move along a smooth curve in three dimensional euclidean 
space — what confines the particle to the vicinity of the curve is a two-dimensional harmonic oscillator 
potential in the plane normal to the curve, of width rj, with its minimum at the position of the curve. It is 
convenient to use an adapted coordinate frame, with coordinates (s, a, /S), where s is arclength along the 
curve, and 77a, r]f3 are distances along the normal n and the binormal b of the curve, respectively. Then 
the position vector r(s, a, (3) of an arbitrary point in the vicinity of the curve is related to the position 
vector R{s) of the curve itself by 

r{s,a,P) = R{s)+-qan{s)+-q[3b{s). (7) 

Taking 77 <C t^^, guarantees that the frame is well defined in the region of physical interest — we do 
assume that k{s) ^ 0. The hamiltonian for the particle, expressed in terms of the adapted coordinates, 
is given by 

H = -—^dAG^''^\dB + V{a,p), (8) 

where A, B range over the adapted coordinates, Gab is the metric induced from the ambient euclidean 
one. Gab = dAr-dsr, G^^ its inverse, and G — {l — rjaK)^ its determinantQ For the confining potential 
we take V{a,(3) = (a^ -t- 0^)/2rf'. Notice that V only depends on the normal coordinates, so that, 
classically, the tangential motion of the particle is free. The normalization condition for the particle's 
wavefunction $ is 

I dsdadl3^/\G\\^\'^ = 1, (9) 

which motivates working with a rescaled wavefunction ^ = jCp^^, obeying 

dsdad(3\^\^ = 1 , (10) 



so that /dad/3|\E'|^ can be interpreted as the probability density for finding the particle at the position 
s along the curve. Accordingly, the hamiltonian undergoes a similarity transformation, 

H^H =\G\^/^H\G\~^I^. (11) 

Plugging in ([5]), and expanding in powers of 77 results in 

H = \h_2 +Ho + 0(7?) , (12) 

V 



where 



Ho = -i(0,-7rL)2-^, (13) 



^We follow here the exposition in I12i — notice that the sign of rj in the expression for G that follows from Eq. (3) in 
that reference is opposite to the one given here — we believe our expression is the correct one. 
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and L ~ i{f3da ~ ad/s) is the generator of rotations in the normal plane. Looking for H eigenstates, 
H'^ = i?^, in the factorized form 

vl/(s,a,/3)=x(«,/3)^(s), (14) 

one is led to consider simultaneous H^2 and L eigenkets Xo" ? 

if_2X<.") = (n + l)xi"\ Lxi")=axi"), (15) 

leading to 

- ^C + iar^/j'^ + i (lar' + cr^r^ - -kA ^j^ = E^^^ (16) 

for the tangential wavefunction (primes denote derivatives with respect to s). Then "^^ — Xa '^<^ ^^^ 
E={n + l)/772 + E,. 

4 Wires With Quantum Memory 

Consider a particle constrained to move on a unit circle, then k{s) — 1 and t(s) — 0. For the normal 
ket, take either of the degenerate doublet |±) = (|10) ±z|01))/\/2, with H-2 eigenvalue 2, and satisfying 
L\±) — ±|±) (we use the standard notation |10) = ajj|00), |01) = al|00), where aj,, oL are creation 
operators of excitations along the axes a, (3). Eq. pB]) then becomes 

(17) 
z o 

with CT = ±1 and solutions 

1 .. ,., zii.2 _ 1 

fcez. (18) 



2/fj' - 1/8 (19) 

(p, (f> denote the standard polar coordinates in the normal plane). 
Consider now a two-parameter deformation of the circle, given by 



- ^c - ^^. 


= Ea^a , 


^ij^) ^ J_e^^s ^ j,i,)__ 


4fc2-l 
8 


esponding 3D states *j_^ = 


= X^^t' 


Xi^pe-^^/^e^'W^, 


4°) = 



3... . .^ , .. .. (20) 



R{s', Ci C) = (cos s, sin s, 0) 

+ C(- cos^ s, sin^ s, 0) + C(0, 0, cos 2s) 

in cartesian coordinates {x, y, z), or 

i?(s;^,C) = (0,-1,0) 

1 (21) 

+ ^(- sin 2s, cos 2s, 0) + C(0, 0, cos 2s) 

in the adapted frame {t,n,b), with ^, C *^ 1 (in the notation of the previous section, ^^ = ^, ^^ = 
C). The corresponding velocity field v — dR/d^\(^=o satisfies the locally arclength preserving condition 
dgV* — Kv" = (similarly for u — 9i?/(9Clc=o)- This enables the physical identification of points of the 
curve with the same s coordinate, for different values of the deformation parameters, so that derivatives 
of the wavefunction with respect to the latter can be meaningfully taken. The above deformation brings 
along changes in k(s) and t{s), which in turn result in the perturbation hamiltonian 

3 

Hi = -C cos 2s + 6<cr (sin 2s ds + cos 2s) , (22) 
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to be added to the zeroth-order hamiltonian in the left hand side of ^T7\ . The resulting first order 
corrected ground state tangential wavefunctions are 

TAi = -^ ("l - (^ ± i3C) COS 2s^ . (23) 



We are still not in the position to use ([4]) though. The reason is that, on the one hand, the true, physical 
wavefunction is the <& that appears in ^, and not the rescaled ^ we have been working with, while, on the 
other hand, our expressions involve the adapted coordinates, which depend implicitly on the parameters 
^ and Cj since the curve being deformed in (PO)) drags with it the adapted frame. Our task then, in 
principle, would be to express the adapted coordinates in terms of the cartesian ones, differentiate them 
with respect to the parameters, treating the cartesian cordinates as constants, and reexpress the result 
in terms of the adapted coordinates. Then the derivatives with respect to the parameters in, e.g., ^, 
would contain contributions both from the exphcit dependence of ^{s,a,/3;^,Q on the parameters, as 
well as the implicit one through the adapted coordinates. Taking the above into consideration, we find 
that the matrix A in ([51) is diagonal, so that the nondegenerate expression for K, Eq. (JH), may be used, 
giving 

i^r--^- (24) 

The fact that the two states in the doublet, corresponding to ct — ±1, pick up opposite geometrical 
phases, allows the detection of the latter by starting out the system in a suitable superposition, e.g., in 
the state |*)t=o = |0, 10) = (|0, +) + |0, -))/\/2, with wavefunction * ~ ae~P^/'^ — a plot of a constant 
probability density surface appears in Fig. [T] Suppose now that the the systen is driven cyclically in the 
^-( plane, tracing the circle ^ = e(l — cos Ai), ( = e sin Ai, with A ^ 1, so that the adiabatic approximation 
is valid. For t = the wire described by (|20|) is a unit circle in the x-y plane, centered at the origin, with 
n radially inwards, and b along z. As t increases, the wire sweeps out the self-intersecting surface shown 
in Fig. [5] (assuming e — .2). The geometrical phase accumulated after each revolution in parameter space 
by the state |0, ±) is ±Aip, where 

A^ = 2ne^Kll^+^=--TTe\ (25) 

so that after m revolutions, the state vector is 

|*)t = (e""'^'^|0,+) +e-""'^"^|0,-))/\/2. (26) 

As a result, the corresponding probability density profile rotates in the plane normal to the curve by an 
angle A0 — Fig. [3] assumes an integer m such that Acf) « tt/2. 

5 Conclusions 

We pointed out a wide class of quantum systems exhibiting geometrical phases: particles constrained to 
move along curves in three dimensional euclidean space, with the shape of the latter in the role of an 
infinite dimensional parameter space. We presented a particular example, involving deformations of a 
circle, and leading to a nonzero Berry curvature, Eq. ([24|) . The effect is observable, through changes in the 
probability density distribution, when superpositions of states are considered, each of which accumulates 
a different, in general, geometrical phase. 

There are obvious directions along which our result might be generalizable: higher dimension and/or 
codimension of the constraint manifold, e.g., the case of surfaces embedded in M.^. Along with such en- 
deavors, a detailed analytical study of the problem we considered here should provide general expressions 
for the curvature in terms of the initial shape of the wire and the velocity fields of the perturbations — 
we defer such an analysis to a longer article, currently in progress. 
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Figure 1: Plot of a constant probability density surface for the state |0, 10), corresponding to excitation of 
the 2D harmonic oscillator along n. The darker colored ring represents the (undeformed) wire. 

A question that we find particularly interesting is the extend to which the details of the deformation 
can be encoded in the probability density distribution. It is clear that, no matter what the initial quantum 
state of the particle is, there are deformations of the wire that pass unnoticed, e.g., those that evolve 
along an arbitrary curve in parameter space and then retrace their evolution along the same curve, back 
to the original shape. This means that one cannot, in general, reconstruct the deformation entirely from 
a comparison of the initial and final probability density distributions, but it is, nevertheless, plausible 
that, in general, more can be inferred than just the fact that there was some deformation. 

Experimental uses of the above effect in quantum wires and nanotubes might be possible, as tiny 
periodic mechanical deformations are encoded cumulatively in a quantum system — this could provide 
sensitive (periodic) motion detectors. In this context, it would be interesting to consider the potential 
amplification of the effect, either in the present, or its higher dimensional variants alluded to above, by 
replacing the single particle employed here by mesoscopic condensates. Another direction worth pursuing 
would be applications to holonomic quantum computing, 

C. C. would like to thank his colleagues at NTUA and, in particular, George Zoupanos and Kon- 
stantinos Anagnostopoulos, for hospitality, and multifaceted support. The authors acknowledge partial 
financial support from DGAPA-UNAM projects IN 121306-3 and IN 108103-3, as well as the EPEAEK 
programme "Pythagoras 11", co-funded by the European Union (75%) and the Hellenic state (25%). 
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